JOURNAL OF SELECTED AREAS IN COMMUNICATIONS, VOL. I, NO. II, NOVEMBER 2008 



1 



Underwater Acoustic Networks: Channel Models 
and Network Coding based Lower Bound to 
Transmission Power for Multicast 

Daniel E. Lucani, Student Member, IEEE, Muriel Medard, Fellow, IEEE, and Milica Stojanovic, Member, IEEE 



00 

o 
o 

<: 
o 

m 



> 

o 
l> 
o 
p 

o 

oo 

o 



X 
S3 



Abstract — The goal of this paper is two-fold. First, to establish 
a tractable model for the underwater acoustic channel useful for 
network optimization in terms of convexity. Second, to propose 
a network coding based lower bound for transmission power in 
underwater acoustic networks, and compare this bound to the 
performance of several network layer schemes. The underwater 
acoustic channel is characterized by a path loss that depends 
strongly on transmission distance and signal frequency. The 
exact relationship among power, transmission band, distance and 
capacity for the Gaussian noise scenario is a complicated one. 
We provide a closed-form approximate model for 1) transmission 
power and 2) optimal frequency band to use, as functions of 
distance and capacity. The model is obtained through numerical 
evaluation of analytical results that take into account physical 
models of acoustic propagation loss and ambient noise. Network 
coding is applied to determine a lower bound to transmission 
power for a multicast scenario, for a variety of multicast data 
rates and transmission distances of interest for practical systems, 
exploiting physical properties of the underwater acoustic channel. 
The results quantify the performance gap in transmission power 
between a variety of routing and network coding schemes and 
the network coding based lower bound. We illustrate results 
numerically for different network scenarios. 

Index Terms — Underwater Acoustic Networks, Network Cod- 
ing, Lower Bound for transmission power, minimal transmission 
power, bandwidth - distance dependence 



I. Introduction 

With recent advances in acoustic communication technol- 
ogy, the interest in study and experimental deployment of 
underwater networks has been growing [1]. However, under- 
water acoustic channels impose many constraints that affect 
the design of wireless networks. They are characterized by a 
path loss that depends on both transmission distance and signal 
frequency, in a far more pronounced way than a terrestrial 
radio system. Thus, not only the transmission power, but also 
the useful bandwidth depend strongly on transmission distance 
[2]. References [6] and [7] present studies of the characteristics 
and design challenges of underwater acoustic networks. 
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In terms of the network layer, a series of routing schemes 
have been proposed for underwater networks over the recent 
years. References [6], [7] and [8] present surveys of different 
routing schemes used in underwater networks. In [3] two dis- 
tributed routing algorithms are introduced for delay-insensitive 
and delay-sensitive applications. Reference [4] presents a 
modification of the dynamic source routing protocol that adds 
location awareness and link quality metrics. In [5] a routing 
protocol based on local depth of the nodes is studied. 

Network coding was introduced by Ahlswede et al [10]. 
Network coding considers the nodes to have a set of functions 
that operate upon received or generated data packets [9]. Work 
in [11] and [12] showed that linear codes over a network 
are sufficient to implement any feasible multicast connection. 
Also, [12] provides an algebraic framework for studying this 
subset of coded networks. Work in [13] presents the idea of 
using linear codes generated randomly in a network. Some 
practical network coding protocols have been presented in [22] 
and [ 14] . Network coding has previously been considered for 
underwater networks, showing better performance than other 
routing schemes. In [15] the problem of error recovery was 
studied in terms of the fraction of delivered packets and total 
number of transmissions. In [16] routing and network coding 
schemes were compared based on the time to complete the 
transmission of a fixed number of packets, and the power 
required to do so. Existing results compare different network 
schemes, but the question remains open as to what is the gap 
between these schemes and the theoretical optimum in terms 
of transmission power 

The objectives of this paper are 1) to establish a tractable 
model for the underwater acoustic channel that will be useful 
for network optimization in terms of convexity, and 2) to 
propose a network coding based lower bound for transmission 
power in underwater acoustic networks by using those models. 
The bound is used to compare the performance of several 
network coding and routing schemes. 

For an underwater acoustic channel both distance between 
two nodes and capacity determine transmission power and 
optimal transmission band. However, the complete model that 
relates these variables is complicated. This paper presents a 
simple closed-form approximation for transmission power and 
optimal operating frequency band as functions of distance and 
capacity. This approximate model stems from an information 
theoretic analysis that takes into account the physics of acous- 
tic propagation, and colored Gaussian ambient noise. 

Reference [2] shows that transmission power as function of 
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distance can be well approximated by P{1) = pP. A similar 
model holds for the bandwidth. The coefficients in this model 
are functions of the required signal to noise ratio (SNR). The 
present work extends this idea of modeling the power and 
bandwidth as functions of distance, but the problem is cast 
into a slightly different framework. Namely, instead of using 
the SNR as a fixed design constraint, link capacity is used 
as the figure of merit. The parameters of the approximate 
model proposed are functions of capacity and distance. This 
approximate model is useful for a broad range of capacities 
and distances. 

The complete model that relates transmission power, trans- 
mission band, distance, and link capacity is provably convex. 
Since the approximate model is used instead of the complete 
model in network optimization, the present work shows the 
operating conditions under which the approximate model is 
convex. 

We assess the minimum transmission power required for an 
underwater acoustic network. A lower bound for transmission 
power is obtained by neglecting interference between the 
nodes and using subgraph selection [9] to establish minimum- 
cost multicast connections with network coding. The convex 
cost function for the network optimization is given by the 
transmission power which depends on the distance and a 
desired data rate via the approximate model for each active 
link. We show that the no-interference assumption in an 
underwater scenario is justified for low multicast rates, and 
randomly placed nodes with inter-node distances less than 
10 km. 

Finally, the network coding based lower bound for transmis- 
sion power is used to compare different routing and network 
coding schemes. We use some of the schemes in [16] for 
a concatenated relay network and extend them to a random 
deployment of nodes in two dimensions. Also, we use an 
ALOHA-like MAC layer instead of a TDMA scheme as in 
[16] since TDMA is not scalable. Furthermore, the problem 
of scheduhng in TDMA is NP-complete [24]. 

The paper is organized as follows. In Section |ll] the com- 
plete model of an underwater channel with Gaussian noise 
is presented, and special characteristics of this channel are 
highlighted. In Section |lll] the approximate model for the 
underwater channel is presented. In Section |IV] convexity of 
the complete model is proven and conditions for convexity 
of the approximate model are studied. In Section [V] the 
lower bound for transmission power using network coding 
and subgraph selection is computed. Section |VI] presents the 
schemes to be used for performance evaluation and gives 
numerical results a two-dimensional network scenario. The 
gap of several routing and network coding schemes to the 
network coding based lower bound is determined. The last 
section summarizes the conclusions of this work and future 
research topics. 

II. Channel Model 

An underwater acoustic channel is characterized by a path 
loss that depends on both distance I in km and signal frequency 
/ as 



Ail,f) = {l/l,,ftaify 



(1) 



where k is the spreading factor, /^e/ is a reference distance, 
and a(/) is the absorption coefficient (Figure 1 in [2]). The 
spreading factor describes the geometry of propagation, e.g. 
k = 2 corresponds to spherical spreading, = 1 to cylindrical 
spreading, and k = 1.5 to practical spreading. The absorption 
coefficient can be expressed in dB/km using Thorp's empirical 
formula for / in kHz [17], which is an strictly increasing 
function. 

The noise in an acoustic channel can be modeled through 
four basic sources: turbulence Nt{f), shipping Ns(f), waves 
Nw{f), and thermal noise Nthif) [2]. The power spectral 
densities (psd) of these noise components in dB re /i Pa per 
Hz as functions of frequency in kHz are presented in [18]. 
These psd's have two important parameters: 1) the shipping 
activity s ranging from to 1, for low and high activity, 
respectively, 2) the wind speed w measured in m/s. Figure 
2 in [2] shows A'^(/) for different values of s and w, and an 
approximation 101ogiV(/) = iVi - r?log(/) for / < 100 kHz, 
where TVi = 50 dB re fj, Pa and 77 = 18 dB/dec. 

Assuming that the channel is Gaussian, its capacity can be 
obtained using the waterfilling principle [25] [26]. The capacity 
of a point-to-point link is 

/ K{l,C) 



df 



(2) 



where B{1, C) is the optimum band of operation and K(l, C) is 
a constant determined to satisfy a given constraint. The band 
B{l,C) could be thought of as a union of non-overlapping 
intervals, B{l,C) = Ui[fl,^i{l,C), f 1^^(1,0')], where each non- 
overlapping interval i has the lower-end frequency flniihC) 
and the higher-end frequency fl^j^Q, C) associated with it. The 
factor 1/ A{1, f)N{f) is shown in Figure 3 in [2] for different 
values of I. This figure suggests that the optimal transmission 
band changes considerably with respect to the distance [2]. 
The transmission power associated with a particular choice of 
{I, C) is given by 



J E 



B{LC) 



Sil,C,f)df 



(3) 



where S{l,C,f) = K{l,C) - A{l,f)N{f)J e B{l,C) is the 
psd of the signal. The corresponding signal-to-noise ratio is 
given by: 



SNR = 



lBii,c)Sil,C,f)A-\l,f)df 
lBii.c)Nif)df 



(4) 



We observe that a choice of (l,C) determines implicitly 
the SNR level. Hence, there is a one-to-one correspondence 
between the pair (l,C) and the pair {I, SNR). The latter 
parameterization was used in [2] to compute the transmission 
power and bandwidth representation to ensure a preset SNR, 
which determines the value of C implicitly. The present 
analysis focuses on using the former parameterization, i.e. on 
determining the power and transmission band that ensure a 
pre-set link capacity. 
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Fig. 1. High and low band edge frequency of the transmission band for 
C = 0.01kbps, k = 1.5, s = 0.5, and w = Om/s. 



Fig. 2. Parameters ai and a2 for P{1, C) and approximate model./ € 
[0, 10] km, C e [0, 2] kbps,fc = 1.5,s = 0.5 and w = m/s. 



A. Dependence on the spreading factor 

The dependence on the spreading factor k is quite simple. 
Let us assume that a model for P{1, C) has been developed for 
a particular value of fc, i.e. P(Z, C, k). To determine P(/, C, k') 
for k' / fc, note that a change in k, the product A{1, f)N(f) = 
{I /Iref)^ o-iff ^ if) constitutes a constant scaling factor with 
respect to /. Therefore, for a link of distance I the term 
B{l,C) will remain unchanged. Thus, if the same capacity 
C is required for k and k' , equation ^ shows that the only 
other term that can vary is K(l,C), i.e. K{l,C,k). Then, 
K{l,C,k') = {l/lr^ff'-''K{l,C,k). Finally, let us use the 
equation (O to determine the relationship between P{l,C,k) 
and P{l,C,k'). 

P{l,C,k')= f (K{l,C,k')~lia{fyN{f))df (5) 
Jb(lc) ^ ' 

= it'^ I (kH, C, k) - itaifyNif)) df (6) 

J B{l,C) ^ ^ 

= li-''P{l,C,k) (7) 

where Im = l/lref to shorten the derivation. Thus, any model 
generated for some parameter k has a simple extension. Also, 
note that the transmission band remains the same for any value 
of k. 

B. Interference Characteristics 

The optimal transmission band of a link was shown to 
change with the distance, under the assumption that the 
channel is Gaussian and that the capacity of a link is obtained 
through waterfilling. If the capacity for a link is low, e.g. less 
than 2 kbps, and the transmission distance is below 10 km, 
the transmission bandwidth will also be low, and its optimal 



location in the spectrum will change dramatically with the 
distance. Figure [T] shows this effect for C — 0.01 kbps. In this 
figure, the high and low band edge frequencies are plotted. 
This figure also shows the high and low band edge frequency 
if an SNR requirement of -20 dB is set, i.e. using the SNR 
instead of the capacity as the fixed parameter As noted before, 
the constraint over the capacity is related to different SNR 
levels depending upon the distance. It is clear that low values 
of C are related to a very low SNR value. 

Figure [T] shows that if two links with the same C — 
0.01 kbps are established, one with / ^ 200 m and the other 
with I ^ 310 m, the optimal transmission bands for these links 
will not overlap; thus, they do not interfere with one another. 
This characteristic of the underwater channel suggests that if 
a network is established in which the nodes are at different 
distances from one another, and each node has a limited range 
of transmission when the data rate requirement is very low 
(all valid assumptions in underwater networks), there will be 
no interference between transmissions of the various links. If 
each link allocates its band optimally, this suggests that a form 
of FDMA is the optimal approach in an underwater network, 
where transmission band is determined by both the distance 
and the required data rate. From a network optimization view 
point, the cost function to be minimized is clearly separable 
under these assumptions, where the channel model for a link 
can be used as the cost function for each of the separable 
terms. 

C. Numerical Evaluation Procedure 

A numerical evaluation procedure similar to that of [2] is 
used to compute the value of P{l,C), B{l,C) and fendihC), 
for a region of values of (/, C). The procedure starts by fixing 
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a target value of the capacity C. Then, for each distance I, 
the initial value of K{l,C) is set to the minimum value of 
the product A{l,f)N{f), i.e. K{l,C) = miuf A{1, f)N{f). The 
frequency at which this occurs, i.e. /o = argmm f A{1, f)N{f), 
is called the optimal frequency. After this, K{1, C) is increased 
iteratively by a small amount, until the target capacity value C 
is met. Finally, this procedure is repeated for each value of C 
in a range of interest. At the n-th step of the procedure, when 
K^^\l, C) is increased by a small amount, the band B^^\l, C) 
is determined for that step. This band is defined as the range of 
frequencies for which the condition A(l,f)N(f) < k'^"-\i,C) 
holds. Then, the capacity c'"-* is numerically determined for 
the current A'(")(;,C) and b("'(Z,C), using the equation 
If c(") < C, a new iteration is performed. Otherwise, the 
procedure stops. 

III. Approximate model 

Evidently, the expressions for the complete model are quite 
complicated to be used in a computational network analysis. 
Also, they provide little insight into the relationship between 
power consumption, B and f^nd 7 in terms of the pair {l,C). 
This motivates the need for an approximate model to represent 
these relations for ranges of C and / that are of interest 
to acoustic communication systems. The model should also 
provide flexibility to changing other parameters, such as the 
spreading factor fc, wind speed w and shipping activity s. 
As shown in Equation (|7]i, any approximate model for the 
transmission power generated for some parameter k has a 
simple extension to any other value of k. Also, a model for 
the transmission band remains the same for any value of k. 

By applying the numerical procedure of the previous section 
for various I and C and fitting the data, it is possible to 



obtain approximate models for power consumption (Eq. |8]i, 
band-edge frequency /e„d(',C) (Eq. [TOl i. and for the band- 
width B{l,C) = fend{l,C) - hm{l,C) (Eq. [HJ. Notc that 
some important properties for these parameters are kept, e.g. 

P{l,Q) = 0. 

P(Z,C) = riC^'lO^D" with (8) 
ai(C) ^ as. + a2C + aiC^ (9) 
a2{C) ^Pi+ P2WlogwC + PiilOlogioiC + l)f 

/e„da,C) = ?°i^'^H0T with (10) 
ai{C) = as + a2mogwC + aiiWlogwCf (11) 
a2(C) = /33 + P2WlogioC + PiiWlogwCf 

B{1, C) = ri^'^ho^T^ with (12) 

ai(C) = Q4 + asWlogioC + a2{WlogioCf + ai{10logiQCf 

(13) 

a2{C) =03+ P2WlogioC + (3i{10logwCf 

The transmission power, band-edge frequency and bandwidth 
of transmission band were computed for a variety of values of 
s, w and two ranges of interest of the pair {I, C): I £ (0, 10] km, 
C G [0, 2] kbps, and I e (0, 100] km, C e [0, 100] kbps. The 
models proposed fitted these cases quite well. Results are 
presented for k = 1.5, w = and s = 0.5, for both cases. 
For the first case, the a and /3 parameters show almost no 
dependence on the shipping activity factor s, especially if the 
wind speed is lu > 0. Thus, the approximate model for this 
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TABLE I 

ai AND 02 APPROXIMATION PARAMETER VALUES FOR P{1, C), fendih C) AND B{1, C), WITH k = 1.5, S = 0.5 AND W = M/S FOR CASE 1: 

I e [0, 10] KM, C e [0, 2] kbps, and Case 2: 1 g [0, 100] km, C e [0, 100] kbps. 



Case 






a2 


Ct3 


ai 


MSE 


01 


02 


/53 


MSE 


1 


P{LC) 


-0.00235 


0.01565 


2.1329 





2.532e-7 


0.014798 


1.0148 


74.175 


5.8979e-5 


1 




4.795e-5 


0.00246 


-0.44149 





3.930e-9 


0.00I7I 


0.07153 


13.738 


3.4706e-5 


1 


B{l,C) 


-5.958e-7 


-2.563e-5 


-0.000305 


-0.30694 


6.599e-9 


-5.I63e-6 


0.33427 


9.6752 


2.9233e-7 


2 


P(l,C) 


-5.617e-5 


0.02855 


2.9305 





0.0001 1 


0.04317 


0.90597 


76.156 


0.000I0I15 


2 




-0.00019 


0.01186 


-0.55076 





I.32e-7 


0.0065157 


-0.032693 


14.739 


7.3024e-5 


2 


S(/,C) 


I.696e-6 


4.252e-5 


-0.00249 


-0.36397 


7.29e-7 


-0.0018252 


0.34788 


10.328 


0.00019414 



case can be simplified to consider w only as part of the model, 
instead of the pair (s,™). 

Figures |2l |3] and |4] show parameters ai and 02 for 
P{l,C), fend{l,C), and B{l,C), respectively, for the first case 
with k = 1.5, s = 0.5 and w = m/s. The values of a's 
and /3's are shown in Table Has Case 1, for parameters ai 
and 02, respectively. These tables also show the mean square 
error (IVISE) of the approximation with respect to the actual 
parameters. A similar result can be found for the second case 
with k = 1.5, s = 0.5 and w = 0. The values of a's and /3's 
are shown as Case 2 in Table [H for parameters ai and 02, 
respectively. For both ranges, the proposed models give a very 
good approximation to the actual numerical values. Also note 
that for the ai(C) parameter of P(l,C), it is possible to use a 
linear approximation, instead of a quadratic model. 

For case 1, the values for a and (3 parameters in the approx- 
imate P{l,C) model show very little dependence with respect 
to s while they show a greater dependency on w. This is not 
unexpected, since the transmission band is at a high frequency 
(between 5 and 40 KHz) where the noise psd is influenced 
more by the w (0{f^^)) than s (0{f^^-^)). Therefore, a further 
approximation is to discard s and consider parameters a and 
(5 to be functions of w only. In particular, a simple model is 
«i = 73 + 72lO/ogio(™ + 1) + -fiiWlogioiw + Vi = 1, 2, 3. 
A similar relation holds for /3i,Vi = 1,2,3. Table HIl shows 7 
parameters for the different a's and /3's. 

IV. Convexity Analysis 

In this section we prove the convexity of the complete model 
of P{1, C) in the entire region of interest, i.e. positive data rates 
and / > 0. Then we discuss the necessary conditions that the 
value of I has to fuUfill to ensure convexity of the approximate 
model. 



A. Convexity of Complete Model 

The convexity of the transmission power of the complete 
model is stated in the following lemma, which is proven in the 
Appendix. Lemma 1 assures that P{l,C) is a convex function 
with respect to C for the ranges of interest of C and / for the 
case of non-overlapping finite bands. 

Lemma 1 P{l,C) is a convex, increasing function with re- 
spect to C, VC > and / > 0, if A{1, f)N{f) > and B{1, C) = 

^^lf^n^{l,C)Jlr,d{^,C)], with < /LdO,^) < C»,V» 

and fUl,C),flUl,C) ^ [fL{l'C),fLa(^,C)],yij^j, i.e. a 
union of non-overlapping finite bands. Proof See Appendix I. 



Another property to be used is given in Lemma 2. This 
lemma assures that if a link between a transmitter i and 
receiver j at a distance I achieves a certain capacity C, another 
node k at distance l' < I from node i, will be able to decode the 
information transmitted from i to j. Note that the transmission 
band is optimal for the link of distance I. 

Lemma 2 C{l,B{l,C)) < C{l',B{l,C)) for l' < I, if 
A{l,f) = {l/lreffcL{f)^ with fc > 1 and a(/) > 1,V/. Proof 
See Appendix II. 



B. Convexity of Approximate Model 

The function P{l,z) represents the minimum power required 
to transmit at a data rate z over a link of distance The 
function P{l,z) was proven to be a convex function with 
respect to z, using / as a parameter (Lemma 1). However, the 
exact model is complicated from a computational viewpoint. 
Let us determine the conditions for which the approximate 
model P{1, z) in equation ^ is convex with respect to z, and 
having I as a fixed parameter We study the case of 2: < 2 
kbps. Since the a and /3 parameters come from fitting the 
data, the only variable left to analyze is the distance /. Note 
that ensuring that P{1, z) is increasing and convex translates 
into the following inequalities: 



ln{l) 



"TTT 



> 



ln{l) 



+ 



in(lO) aa2(z) 



+ 



(14) 



(15) 



There is both a linear and a quadratic constraint upon / to 
ensure convexity. Since these constraints are also functions of 
2, the range of values of this parameter should be considered. 
From previous results for the fitting parameters, it is possible 
to determine some properties of the model for 2 < 2 kbps. In 
terms of the parameters of interest, ai < 0, a2 > 0, 2aiC -f 
a2 > 0, /3i > and (32 > 0. Thus, for the choices of ai (2) and 
02(2), the first and second derivatives of these functions with 
respect to 2 are di{z) > 0, 01(2) < 0, 02(2) > 02(2) < 0. 
Using these conditions, the constraints (fl4] i and (fTSl l can be 
simplified to 



ai(z) 



2di{z) 



in(io) / d2{z)d\{z) _ 02(2) 



di{z) 



di{z) 



+ 



d'ljz) 
4ai(2)^ 



(16) 
(17) 
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TABLE II 

Approximation PARAMETERS of a and /3 for P(/, C), I e [0, 10] km, 
C G [0, 2] KBPS, k = 1.5,s = 0.5 





71 


72 


73 




5.2669e-6 


-0.000157 


-0.004575 




-2.97Ie-5 


0.000865 


0.029306 


013 


0.000152 


0.01809 


2.4586 


Pi 


9.924e-6 


-0.00027 


0.012288 


(il 


7.799e-6 


-0.000219 


1.0118 


/33 


0.068091 


1.3659 


73.144 



where the term under the square root is positive which ensures 
real values of I. Note that for each value of z there is a 
minimum value of /. Let us use the values of Case 1 in Table U 
to determine the (l,z) region for which the approximate model 
is convex. For these values, if the distance between to nodes I 
is at least 13 m, for any value of z <2 kbps the model will be 
convex. The limitation to I >13 m is related to the sampling 
of the distance used for computing the parameters of the 
approximate model. For all practical purposes the approximate 
model P{l,z) is convex. 

V. Lower Bound to Transmission Power in 
Underwater Networks 

The problem of achieving minimum-energy multicast us- 
ing network coding in a wireless network has been studied 
previously [9]. A wireless network, as presented in [9] can 
be represented through a directed hypergraph H = C^jA) 
where K is the set of nodes and A is the set of hyperarcs. 
A hypergraph is a generalization of a graph, where there are 
hyperarcs instead of arcs. A hyperarc is a pair {i, J), where i, 
the start node, is an element of K ,and J is the set of end nodes 
is a nonempty subset of A. Each hyperarc {i, J) represents 
a broadcast link from node i to nodes in the nonempty set 
J. Let us denote by Zij the rate at which coded packets are 
injected into hyperarc (i, J). If the cost function is separable, 
the optimization problem can be expressed as follows 



min ^fi^^J/^) 
{i,J)eA 

subject to z e Z 

„(*) 



{j\{i,j)eA}jeJ 
(*) 



E 

{3\U-i)eA},i£i 



(t) 



xl>^>Q,\/{i,J)£A,j^J,t£T 



with 



r(t) 



R if j = s, 
"R if j = t, 
otherwise 



(18) 



(19) 



where T is a non-empty set of sink terminals, a source s, a 
multicast rate 7?, and a fixed transmission duty cycle at each 



link 9. xfj- represents the flow associated with terminal t, sent 



through hyperarc (i, J) and received by node j e J. 



In the underwater scenario this formulation is used to 
establish a lower bound on the transmission power required 
to achieve a multicast rate R. Assuming no interference for 
transmissions in different hyperarcs yields a separable cost 
function. Note that if interference was taken into account, the 
power to reach the desired data rate would increase. Then, the 
cost function f{zij) for each particular hyperarc corresponds 
to a link transmission power P{l,Zij) in order to obtain the 
minimum transmission power required to achieve a data rate 
of Zij, where I represents the distance from i to the farthest 
node j e J. For the lower bound computation, continuous 
transmission (6* = 1) is assumed. A simplification of this 
problem can be made under the assumption that transmissions 
are omnidirectional, and considering the fact that if a node 
transmits over a certain range, all nodes in that range will be 
able to receive the information. This was proven in Lemma 
2. Finally, the model for this channel ensures that any value 
of Zij can be achieved if enough power is used. Thus, the 
constraint set Z can be dropped. 

Although the problem for minimum-cost multicast is well 
known for wireless radio networks, the cost function pre- 
sented here is different because it represents the minimum 
transmission power for an hyperarc transmitting at a data 
rate Z, which is given by the power needed to transmit at 
capacity C = Z, without assumption on technology or, more 
importantly, a specific transmission band which is usually the 
case for wkeless radio networks. Thus, we are providing a 
lower bound valid for any acoustic underwater network for 
the case of Gaussian noise. 

VI. Performance Comparison 

For this study, five schemes are considered. The first scheme 
corresponds to the lower bound to the transmission power 
using network coding given by solving the problem in Section 
|V]with 6 = 1. The second scheme corresponds to solving the 
problem in Section |V] for 6* < 1, in order to study the effect 
of using a duty cycle for link transmissions in underwater 
networks over interference and transmission power. The third 
scheme corresponds to using the paths chosen by the optimal 
scheme but establishing a SNR requirement for the transmis- 
sion links with the objective of studying interference when 
the SNR requirement is increased. The links are considered 
to transmit continuously. The schemes (4) and (5) consider 
implementations of network coding in a rateless fashion with 
the implicit acknowledgment (ACK) [16] and routing with 
link-by-link ACK using an ALOHA-like MAC layer Let us 
explain in more detail each of the schemes. 

l)Network coding based lower bound to transmission 
power: Transmission power is computed by solving the con- 
vex optimization problem in |V] and it provides a lower bound 
on the optimal transmission power for networks operating at 
low data rates. For this computation, continuous transmission, 
i.e. a duty cycle of 6* = 1 is used. This scheme is used as the 
gold standard to which the remaining schemes are compared. 
The no-interference assumption is assessed by computing the 
average percent of the randomly deployed networks that have 
at least a link which suffers from severe interference, i.e. a 
signal-to-interference ratio (SIR) below 3 dB. 
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Fig. 6. Subgraph selection for scheme (4). Selected subgraph with val- 
ues of Zij to provide a unicast rate of 100, where dashed lines repre- 
sent unused transmission ranges. Note that hyperarcs include in the same 
range, e.g. the possible hyperarcs with node 1 as the staring point are 
1{2},1{2,4},1{2,3,4}. 



2) Network Coding with optimal power consumption for 
links with fixed duty cycle: Transmission power is computed 
by solving the convex optimization problem in (|V]i for links 
with a fixed duty cycle, i.e. 9 < 1. This value provides a 
lower bound on the optimum power consumption for networks 
operating at low data rates when links have a particular duty 
cycle. By convexity of the cost functions used, this bound 
will be higher than for the previous scheme. This scheme 
is used to illustrate the effect upon transmission power and 
interference when the links transmit at a fixed duty cycle 6 < 1 
by comparing this scheme to the previous scheme. 

3) Network Coding with SNR requirement on link trans- 
mission: This scheme is a heuristic scheme that imposes 
an SNR requirement for transmissions. Using the subgraph 
selected by solving the problem in (|V]i for 6 = 1, it computes 
the SIR on the different links for a variety of SNR constraints 
using the models of transmission power and band in [2]. As 
in scheme (1) and (2), the percent of randomly deployed 
networks with at least a link with severe interference is 
computed. Results of this scheme suggest that continuous 
transmission with a moderate SNR requirement causes severe 
interference. A solution to this problem is to use of a duty 
cycle 6* < 1, similarly as in scheme (2) when there is an SNR 
requirement. 

4) Network coding in rateless fashion with implicit ACK: 

For a concatenated relay network as in Figure |5] the path 
between a source node and sink node is fixed and includes 



all relay nodes. If a node b is closer than node a to the 
collecting node, a is said to be upstream with respect to 
b, and b is said to be downstream with respect to a. For 
the concatenated relay network this ordering is quite natural. 
This problem was studied in [16]. For a two-dimensional 
scenario, subgraph selection [9] with linear and separable 
cost functions are used to determine the active Unks in the 
network and the transmission power required for each link. 
The cost function of each hyperarc is computed based on the 
approximate formulas for transmission power and bandwidth 
for a fixed SNR level [2]. The weight of each link is given 
by DP{l,SNR)/B{l,SNR), where D is a constant common 
to all Unks related to the number of transmitted bits per 
burst and modulation used. For the performance computation 
using optimal modulation, i.e. Gaussian signaling, the weight 
for each link in the path is DP{1, SNR)/C{1, SNR), where 
C{l,SNR) is the function of capacity related to the pair 
{I, SNR). We assume that the coding is over very large number 
of data packets. This could be extended computing an error 
probability based on error exponents [26]. Once the subgraph 
has been selected, if several links share the same transmitting 
node, this node will randomly choose the link to use. The 
weight of each link in the random choice is given by the 
fraction of data rate the optimization problem assigned to each 
of these links. Let us consider the network in Figure |6] as an 
example, where dashed lines represent unused hyperarcs. If 
node 1 transmits to both node 2 and 4, but we send a rate of 
90 units through hyperarc 212, while we send 10 units through 
hyperarc 212,4, then when node 1 transmits it will do so 90 % 
of the time to reach node 2 only, and 10 % of time using 
enough power to reach nodes 2 and 4. Finally, we have to 
determine which nodes are upstream and downstream to each 
node in the subgraph. If the subgraph corresponds to a single 
path, the choice is clear. If there are multiple paths, we use the 
following heuristics: we start by ordering the nodes starting at 
the transmitter and looking at the nodes directly connected to 
it in the optimal subgraph. These nodes are ordered as follows: 
the node associated with the link with higher data rate from 
the transmitter is considered to be directly downstream from 
the source node, the node with the second highest data rate is 
considered to be downstream with respect to the previous one, 
and so on. In Figure |6l 2 is directly downstream of 1, and 4 is 
downstream of 2. Once all nodes connected to the transmitter 
(let us call this set of nodes S) are ordered, we proceed to order 
the nodes connected to 5 by a similar procedure as for the case 
of one node described before. In the example, S = {2, 4} and 
the nodes connected to it are {3, 4}. If a node connected to 
one of the nodes in S has already been ordered, like node 4 in 
the example, the link is discarded keeping the previous order 
of the nodes. We update S with the nodes that were connected 
to S and not previously in it, until we reach the receiver. For 
the network example in Figure |6] the ordering is 1,2,4,3. 

For this particular scheme, once a relay node gets its 
first coded packet, i.e. a packet formed by a random linear 
combination of data packets, it will transmit until the receiving 
node sends a confirmation that all the information has been 
received. The same happens at the source node. However, 
nodes eavesdrop on other transmissions. If a node receives a 
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coded packet from a node further downstream with the same, 
or a greater number of degrees of freedom than what it has, it 
will stop transmitting and update its information if necessary. 
Degrees of freedom in this setting represents the number of 
packets that were linearly combined to form the coded packet 
as in [16]. The node will resume transmitting if an innovative 
packet, i.e. a packet with a new random linear combination of 
data packets useful for decoding the information, is received 
from a node upstream. The sink node will retransmit a coded 
packet with its degrees of freedom when a coded packet is 
received. This strategy assumes that there is a mechanism that 
informs the collecting node about the number of degrees of 
freedom that constitute the total message or that this number 
is fixed a priori. 

5)Routing using link-by-link acknowledgement: For a 
concatenated relay network, the path between the source node 
and the sink node is fixed and includes all relay nodes. This 
problem was studied in [16]. For a two dimensional scenario, 
the sink and the source are chosen randomly and the shortest 
path is computed before starting data transmission in unicast. 
The weight of each link is computed based on the approximate 
formulas for transmission power and bandwidth for a fixed 
SNR level in the same fashion as the cost function per link 
of scheme (4). In the current scheme, every time a node 
receives a packet, it will retransmit the packet and send an 
acknowledgement to the previous node. Once a packet has 
been acknowledged, the node can start transmitting a new data 
packet in its queue. If it has no new packets to transmit, it 
will only transmit if a node upstream sends new information, 
or sends a previous packet, in which case the node will 
acknowledge this packet. 

In terms of the physical layer, schemes (4) and (5) use 
both PSK modulation, which implies the use of a data rate in 
each link that is lower than capacity, and Gaussian signaling 
assuming that the encoding is over a large number of bits. 
In order to deal with the SNR requirement, we use an 
approximate model for the transmission power, high band edge 
frequency and bandwidth as functions of SNR similar to the 
work in [19]. When PSK modulation is used, probability of 
packet error due to noise over the link from node i to j is 
obtained from the probability of bit error by Ppacket Eitoi (*i j ) = 
1 — (1 — Pbit eiTor)"> where n is the number of bits in the packet, 
and Pbit error IS computcd using the standard PSK bit error 
probability. Note that nodes farther away from the transmitter 
have some probability of receiving the packet correctly. For 
Gaussian signaling, the probability of packet error due to noise 
is considered to be zero for all nodes in range, and 1 for all 
nodes further away. 

In terms of the MAC layer, schemes (4) and (5) use an 
ALOHA-like MAC layer. This ALOHA protocol considers a 
fixed number of bits per data packet and uses the optimal 
transmission band for an SNR requirement per link. Thus, the 
duration of the transmitted packet depends on the transmission 
distance [2]. Every node has a probability to access the 
medium every T units of time following a Bernoulli process. 
Transmission delay is considered using a typical value of 
sound speed (1500 m/s). Figure |5] shows an example of using 
this MAC layer for three nodes with Di >> D2- In this 
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Fig. 7. Percent of deployments in a fixed square of 5x5 km^ witli SIR < 
3 dB in at least one link vs number of nodes deployed in that area, for the first 
three schemes. For scheme 19 = 1, while scheme 2 is shown with different 
values of 8 to achieve unicast rate of R = 0.1 kbps. Performance for scheme 
3 is shown for different SNR values. 



example, when node 1 transmit a packet to node 2, this 
packet also reaches node 3. Note that the duration of the 
packet transmitted from node 1 to node 2 (Packet A) is large 
compared to the packet transmitted from node 2 to node 3 
because of the relation of distance to bandwidth/capacity for 
a fixed SNR value mentioned above [2]. Once node 1 has 
transmitted the packet it will try to transmit again, and it has 
some probability to start transmission every time slot T. Let 
us assume that node 2 has a data packet for node 3. Figure 
|5] shows the case when the data packet transmitted from node 
2 to node 3 suffers a collision at node 3 with a new packet 
transmitted from node 1 to node 2. We consider that a collision 
at any receiver causes a loss of all packets involved in the 
collision for that receiver 

Let us study some numerical results that correspond to a 
network in which nodes are deployed randomly in a two 
dimensional space. Unicast connections of rate 7? are es- 
tablished, i.e. the network has one transmitter, one receiver 
chosen randomly, and, possibly, several relay nodes. The 
number of nodes ranges from 3 to 8. The transmission power 
lower bound, as an average over random deployments, will 
be compared with transmission power of schemes for routing 
and network coding. Also, a comparison between the schemes 
(1) and (2) in terms of interference is presented. Note that the 
transmission power lower bound is computed assuming that 
all nodes are within the transmission range of the others, i.e. 
full connectivity. 

Figure |7] illustrates the effect of introducing a duty cycle 
6 (dashed lines) for link transmission under a random de- 
ployment in a 5 X 5 km^ square. For 6 = 1 in Figure |7] 
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Fig. 8. Average transmission power of networks deployed randomly on a 
1x1 km^ square. Network schemes operating at SNR = 10 dB and a lower 
bound for transmission power (scheme 1) is presented. Model used considers 
k = 1.5, s = 0.5 and ui = m/s. 
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Fig. 9. Average transmission power of networks deployed randomly on a 
1x1 km^ square. Network schemes operating at SNR = 10 dB using Gaussian 
signaling. Model used considers k = 1.5, s = 0.5 and w = m/s. 



which corresponds to scheme 1, note that less than 3 % of the 
random deployments cause severe interference (SIR < 3 dB) 
over at least one link. This corroborates the no-interference 
assumption used during the analysis to obtain a lower bound 
for transmission power Furthermore, this percentage seems 
to have little dependence on the number of nodes deployed 
in the network. When a value of 6* < 1 is used. Figure |7] 
shows that the percentage of deployments with SIR < 3 dB 
increases for the initial decrements of 6, but decreases as 6 
becomes very small (below 1.5% for 9 = 0.01). Although this 
may seem counter intuitive, introducing a duty cycle causes 
the link to transmit at a higher data rate when it is active which 
translates to using more bandwidth and power to achieve 
that data rate in the underwater channel. Although duty cycle 
reduces interference by not using the channel continuously, the 
combined effect with the increased transmission bandwidth 
and power causes more interference for initial decrements 
on the value of 6. This is a transient effect, and it has a 
breaking point for a small value of 6 when the probability 
of having interference is small. As the value of 9 decreases 
the transmission power can be shown to increase. This is an 
expected effect since the cost function is convex and the value 
of 61 is a constant parameter to all links in this problem. 

Figure |7] shows the results in continuous lines for the third 
scheme with different SNR requirements. The figure presents 
the percentage of random deployments that have at least one 
link suffering from severe interference. For very low SNR, 
the assumption of no-interference is justified. However, even 
for SNR = -5 dB the percentage of deployments with severe 
interference for a unicast connection increases dramatically, 
especially when the number of nodes in the network increases. 



A similar effect occurs when SNR = dB. One way to reduce 
interference while having an SNR requirement is to use a 
similar approach as is scheme (2), i.e. to have a transmission 
duty cycle in each of the links. Schemes (4) and (5) show 
an implementation using an ALOHA MAC protocol, where 
every link has an associated duty cycle when it has some data 
to transmit. 

Let us compare the transmission power of scheme (4) and 
(5) to the lower bound using both PSK and Gaussian signaling 
in a 1 X 1 km^. Figure [8] shows the average transmission power 
for different number of nodes in the network, both active and 
inactive, i.e. before determining the shortest path or solving 
the subgraph selection problem, with a transmission power 
computed to obtain a burst SNR = 10 dB. The average data 
rate for the different schemes is, R ^ 1 kbps. This figure 
shows optimal signaling (Gaussian signalling) and a PSK 
modulation, which illustrates that close to 6 dB in the gap 
between a PSK modulation and the lower bound is due to 
the choice of the modulation. Notice that the gap between the 
average transmission power for Gaussian signaling and the 
lower bound of i? = 1 kbps in Figure [8] is about 1 1 dB for 
scheme (4) and 13 dB for scheme (5). Also, it shows that 
this gap is maintained as more nodes are deployed. Some part 
of the gap is related to the MAC protocol used. Another is 
related to the 10 dB SNR requirement which is usually used 
for a practical implementation. 

Figure |9] compares transmission power for different data 
rates using Gaussian signaling. The number of transmitted 
bits was kept constant, while the transmission probability over 
each link was increased to achieve the desired rate. Note that 
transmission power increases by 3 dB for scheme (4) while it 
increases by almost 5 dB for scheme (5) when the data rate is 
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Fig. 10. Energy for transmission for networks deployed randomly on a 
1x1 km^ square. Network schemes operating at SNR = 10 dB using Gaussian 
signaling. Model used considers k = 1.5, s = 0.5 and w = m/s. 



increased from 1 kbps to 2 kbps, i.e. the gap between scheme 
(4) and (5) increases as data rate increases. This figure shows 
also that the gap between schemes (4) and (5) is very low 
when the data rate is 0.2 kbps. Note that an increase in data 
rate is related to an increase in the collision probability in 
the ALOHA protocol. For the same setting. Figure [TO] shows 
transmission energy for both schemes. The energy required 
for transmitting at the chosen data rates remains constant in 
the case of network coding, while it increases for routing 
when high data rates are attempted. Multiple transmissions 
of one packet are the main cause of the increased energy 
consumption for scheme (5), caused both by packet losses 
and long delays in transmitting an ACK packet given the 
ALOHA MAC layer. While scheme (4) transmits innovative 
packets at each transmission, scheme (5) tries to retransmit 
the same packet if no ACK has been received. Consider the 
case of a long delay in transmitting an ACK, i.e. the packet 
was correctly received but the ACK is transmitted a long time 
after reception, scheme (5) can generate several transmissions 
of the same data packet. This involves an additional energy 
consumption. For the same number of transmissions, scheme 
(4) will transmit several innovative packets, which are useful 
in decoding the information at the receiver 

These results illustrate that coding, subgraph selection and 
the eavesdropping capabilities associated with network coding 
allow a better performance when the collision probability 
increases. However, notice that when transmission rates are 
low the benefits of network coding are less marked. This is 
explained by the fact that an implicit ACK might or might 
not be received by an upstream node. If it is not received, the 
node will keep transmiting innovative packets. This effect is 
particularly evident when Gaussian signaling since only nodes 



in range of transmission will correctly receive a packet. If we 
use a similar example as in figure |5] node 1 will not receive 
any implicit ACK from node 2, and it will continue to transmit 
until informed that all information was received at the sink 
node. Thus, an explicit ACK procedure should be used if node 
deployments are likely to produce these situations, especially 
if only one node is actively generating new data packets. 

VII. Conclusions 

A tractable model for the underwater acoustic channel has 
been established and used for network optimization. A lower 
bound for transmission power in underwater acoustic networks 
was obtained based on network coding, and this bound was 
used to compare the performance of several network coding 
and routing schemes. 

The closed-form approximate models for the time-invariant 
acoustic channel where shown to provide a good fit to the 
actual empirical values by numerical evaluation for different 
ranges of distance / and capacity C, as well as noise pro- 
files corresponding to different shipping activity factor and 
wind speed. The parameters obtained for these approximate 
models can be used in the case that a different spreading 
factor is needed, since the band-edge frequency fendi^C) 
and the bandwidth B{l,C) were shown to be invariant to the 
spreading factor k, while the power scales as P{l,C,k') = 
{l/lref)^ ~^P{h C, k), where / is in km. Also, the approximate 
model of P{l,C) was shown to be almost independent of the 
shipping activity factor s while having a marked dependency 
on the wind speed w for I < 10 km and C < 2 kbps. This 
dependence on w is quite smooth and can be approximated by 
a simple model, resulting in a complete model for the P{1, C) 
for a range of values [l,C) that is of interest to a typical 
underwater communication system. 

We show that the complete model is convex for all C > 
and distances I > 0. Since the complete model is compli- 
cated to solve a network optimization problem, we present 
conditions on the minimum distance between nodes that ensure 
convexity of the approximate model. Convexity of this model 
allows us to use it in more complex scenarios, for example, 
in the framework of layering as optimization decomposition 
[20] [21]. 

This work shows that the no-interference assumption used in 
the computation of the lower bound to transmission power in 
the underwater scenario is justified for low multicast rates, and 
randomly placed nodes with distances under 10 km between 
each other We present numerical results that confirm the 
validity of this assumption, showing that less than 3% of the 
links suffer severe interference. We show that solving the op- 
timization problem when the links have a fixed duty cycle for 
transmissions can reduce interference if the duty cycle is low 
enough. However, high duty cycles can effectively increase 
interference because of the high dependence of the bandwidth 
to the transmitted data rate in an underwater channel. 

The network coding based lower bound was used to de- 
termine the gap of different medium access protocols and 
network schemes for some multicast rate in underwater net- 
works. This comparison is carried out for several routing 
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and network coding schemes. We present numerical results 
in unicast scenarios by deploying nodes in a two dimensional 
environment. 

Network coding with implicit acknowledgements introduced 
in [16] has better performance than routing with link-by- 
Unk ACK in terms of transmission power, especially when 
the probability of collision is increased. The gap between 
routing with hnk-by-hnk ACK and network coding with 
impUcit ACK in terms of transmission power was shown to 
increase as the transmission probability increased, which is 
closely related to the probability of packet collision. However, 
study of a network coding scheme with an expUcit ACK 
should be considered to improve performance. Also, network 
coding with impUcit acknowledgement compared to common 
rateless network coding schemes allows to save resources, e.g. 
memory required in the nodes, and rate adaptation following 
a similar analysis as in [23]. Also, there has been previous 
work on network coding based Ad-Hoc protocols, such as 
CODECAST [22], which could be extended to use implicit 
acknowledgements, and adapted to the underwater acoustic 
channel. 

Future research should consider using error exponents [26] 
to compute the error probability when Gaussian signaling is 
used. This provides a more accurate estimate of transmission 
power when encoding is performed over a finite number of bits 
with a data rate approaching capacity. Finally, future research 
should address the issue of scalability, i.e. study if the results 
are valid when the number of nodes in the network is increased 
but maintaining the same node density. 

Appendix I 
Proof of Lemma 1 

We consider a set H of bands, each band 
i e H having a flndi^^^) flnii^^C!) associated 

to it. Then, P{l,C) = Y.iK{l,C)Uindil,C) - 
fLil,C)) - T,i!^£^^c) A{lJ)N{m and C = 

4'"(^c? ^""3^ {-WlWUl) Using the Leibniz Integral 
rule, ffle fact that K{l,C) = ^(^ /e„d(^, C'))iV(/e„d(Z, C)) and 
K{l,C) = A{l,f,„,{l,C))N{f,„,{l,C)), that A{l,f)N{f) is 
independent of C, and that the derivative of the sum is the sum 
of the derivatives, =2I^(i^Y,ifind{l,C)-fLdl,C)). 

i 

Taking the second derivative: 

= H '^^UinAl,C) - fLi{l,C)) (20) 

i 

4-T ( 'lI^ndhEl 2Jln££l\ (2)) 

^ l^i — I 5C 5C ) ^'^^> 

Taking the derivative of C with respect to itself, using Leibniz 
Integration Rule and K{\,C) = A{lJ^^a{l,C))N{f,^S,C)) 
and K{l,C) = A{l,h„i{l,C))N{fi„i{l,C)), then: 

1 = /n(2)A(/.6-) E {uLdil,C) - fUUC))) {IT) 

i 

Since K(l, C) > for any I > and C > by the physics of 
the channel and fi^^{l, C) - fl^iil, C) > 0, VC > 0, / > , and 
the i bands are non/overlapping and ln{2) > this implies that 



^^^4^ > 0. Then = ln(2)K(l,C) > 0,V/ > 0,C > 0. 

Taking a second derivative to the C expression with respect 
to itself: 

{^^yj:^ifindii,c)-n„s,c))= (23) 

Y.^flnd(h C) - fUl, C)) (24) 

i 

Thus, = {^)'Y.iilnai.l,C)-iUl,C)) where 

(/end(''C)-/L(^C')) > 0,VC > 0, finite and non-overlapping 
and^i^>O.Thus,^^^>0 

Appendix II 
Proof of Lemma 2 

Since A{IJ) = {l/lrei?a{ff and = 
{k/lref){l/lreff-^a{f)^ + {I /lreffln{a{f))a{ff > sincc 
a{f) > 1 and l„.y > 0. Then, A{l,f) > A(l',f),l > I'. 
Also, A{if)Nlf) ^ 1'^/ £ B{l,C) which implies 
^og2( ^(fj)'^(j) ) > 0,V/ € B{l,C). Let us compute the 
capacity of a link of distance l' when we use the optimum band 
and spectral density for a link of distance / and capacity C, 
i.e. B{l,C) and S{l,C,f) = K{l,C) - A{1, f)N{f), f € B{l,C), 
respectively. Then, 

C{l',B(l,C)) = /^(,^) log, (l d/(26) 

> lBii,c) 1092 (l +ML^)^-1^) df = C (27) 
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